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摘 要 : 我 们 利用 Lomb-Scargle periodogram 方法 和 加 权 小 波 Z 变 换 (Weight Wavelet 
Z-transform) 7 & 47 T #2 EKCGRaBS J0929+5013 42 了 J2146-1525 在 15 GHz 射电 和 伽 马 波段 
约 12 年 的 观测 数据 。 分 析 结 果 表 明 炮 变 体 CGRaBS J0929+5013 和 J2146-1525 在 射电 波段 可 能 
存在 准 周期 振荡 信号 。 它 们 的 准 周期 分 别 为 1435 天 (3.93 年 ) 和 1321 天 (3.62 年 )， 且 置信 度 
均 >3O (99.7%); 而 在 伽 马 波 段 光 变 周 期 性 较 弱 。 基 于 超大 质量 双 黑 洞 系统 中 的 螺旋 喷 流 模 


Al fe T 42 & AACGRABS J0929+5013 和 J2146-1525 的 主 黑洞 质量 M 的 分 别 为 4.3x10?47。 和 


2.7x10°M 。。 为 了 探讨 射电 和 伽 马 波段 的 光 变 起 源 ， 我 们 用 离散 相关 函数 分 别 计算 了 雁 变 


体 CGRaBS J0929+5013 和 了 2146-1525 射 电 与 伽 马 波 段 之 间 的 相关 性 ， 发 现 它们 之 间 没 有 明显 
的 相关 性 。 这 表明 它们 的 辐射 可 能 来 自 不 同 的 区 域 。 
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活动 星系 核 (Active Galactic Nuclei, AGN) 是 一 类 特殊 的 星系 ， 其 上 存在 超过 恒星 内 部 核 
反应 的 产能 、 非 热 辐 射 、 高 能 X 和 伽 马 射线 、 物 质 的 喷射 等 猛烈 的 活动 现象 和 剧烈 的 物理 
过 程 , HPR Blazar) 是 活动 星系 核 的 一 个 子 类 。 HEE KM AMR BL 型 天 体 (BL 
Lac) 和 平 谱 射电 类 星体 (FSRQ) ，BL Lac 天 体 和 FSRQ 天 体 的 区 分 依据 是 发 射线 的 强 弱 ， 
发 射线 比较 弱 或 没有 的 是 BL Lac 天 体 ， 比 较 强 的 是 FSRQ APR. WEARER PE AD 
变 ， 比 较 高 的 辐射 偏振 和 非 热 连续 光谱 等 特征 中 。 绝 大 部 分 的 泡 变 体 的 光 变 具有 随机 性 ， 
而 具有 准 周期 性 光 变 的 次 变 体 极 少 , 但 是 也 有 少数 准 周期 候选 体 被 文献 所 报道 。 随 着 观测 技 
术 的 提升 和 数据 的 积累 ， 不 断 有 新 的 候选 体 被 发 现 。 其 中 最 著名 的 是 Sillanpaa 等 人 1988 
年 发 现 OF 287 FENG V 波段 具有 大 约 12 年 的 周期 性 爆发 四。 这 一 周期 性 爆发 现象 可 以 用 
超大 质量 双 黑 洞 模型 很 好 地 解释 。 四 年 后 ，Kidger M 等 人 BI 分 析 了 OJ 287 光学 B 波段 一 百 
多 年 的 数据 ,发 现 B 波段 爆发 周期 为 11.6 土 0.5 年 , 证 实 了 Sillanpaa 等 人 中 的 结论 ， 并 预测 
了 下 一 次 爆发 最 大 值 出 现在 1994 年 8 月 左右 。 随 后 越 来 越 多 的 科研 工作 者 参与 准 周期 搜寻 ， 
Ackermann M 等 人 四 发 现 PG 1553+113 在 光学 和 伽 马 波段 具有 2.18 +0.08 年 的 准 周期 震荡 ， 
Penil P 等 人 GI 用 多 种 方法 分 析 Fermi Large Area Telescope (LAT) 数 据 , AH 11 个 活动 星系 核 
存在 准 周 期 震荡 ， 其 中 9 个 是 首次 发 现 ，Yang Xing 等 人 器 发 现 3FGL J0449.4-4350 在 光学 
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和 爷 马 波段 存在 450 天 左右 的 准 周期 震荡 , 基于 超大 质量 双 黑 洞 模型 的 基础 上 估算 了 其 主 黑 
洞 质量 ， 约 为 7.7x10?M7。， 约 云 露 等 人 吕 发 现 汐 变 体 CGRaBS J0835+6835 的 射电 波段 存在 


560 天 左右 的 光 变 准 周期 ， 并 用 螺旋 喷 流 模型 估算 了 主 黑 洞 质量 在 3.7X 10s M 。~3.5X10? 


MM。 之 间 。 


在 欧文 斯 谷 射电 天 文 台 (Owens Valley Radio Observatory , OVRO) 40m 射电 望远镜 和 
费 米 伽 马 射线 太空 望远镜 (Fermi Gamma-ray Space Telescope) 的 观测 数据 库 中 , 我 们 用 时 间 序 
列 分 析 法 进行 准 周期 信号 的 搜索 。 搜 索 结果 表明 BL Lac 天 体 CGRaBS J0929+5013(4FGL 
J0929.3 +5014) 和 FSRQ 天 体 CGRaBS J2146-1525(4FGL J2146.4-1528) 可 能 存在 准 周期 震荡 ， 
因此 ， 我 们 对 J0929+5013 和 J2146-1525 射电 和 伽 马 波段 的 数据 做 了 进一步 分 析 ， 探 究 两 个 
波段 之 间 的 光 变 相关 性 以 及 辐射 机 制 。J0929+5013(a=09 29 15.4400532384,5=+50 13 
35.987058917) 是 一 颗 z( 红 移 )=0.37039 的 天 体 ， 在 B、V 波段 的 流量 分 别 为 0.00049 mJy, 
0.00069 mJy], J2146-1525(a=21 46 22.97933395,5=-15 25 43.8856129) 是 一 颗 z( 红 移 )=0.70 中 
的 天 体 ， 在 B、V 波段 的 流量 分 别 为 0.00030 mJy, 0.00045 mJyls1。 本 文 的 结构 如 下 : 第 一 
章 介绍 研究 所 用 的 数据 ; 第 二 章 介 绍 Lomb-Scargle periodogram (LSP)! IF IE AM GR 
第 三 章 介 绍 加 权 小 波 Z 变换 (Weighted Wavelet Z Transform, WWZT)0519 方 法 和 检测 结果 ; 第 
四 章 分 析 了 两 个 波段 的 相关 性 ; 第 五 章 为 结论 与 讨论 。 


1 数据 分 析 


我 们 分 别 从 欧文 斯 谷 射电 天 文 台 (Owens Valley Radio Observatory, OVRO) 40m 射电 望 远 
镜 的 官网 (Chttps:/sites.astro.caltech.edu/ovroblazars) 和 费 米 伽 马 射 线 太 空 望 远 镜 官网 
(Fermi-LAT, https://fermi.gsfc.nasa.gov/ssc/data/access/lat/10yr_catalog/ap_lIcs.php) HY XW Wl) Be He 
中 收集 了 JO929+5013(J0929.3+5014) All J2146-1525(J2146.4-1528)15 GHz 射电 波段 和 0.1-200 
GeV MERKRA 12 年 左右 的 数据 .OVRO 40 m 望远镜 配备 了 离 轴 双 光束 Dicke-switched A 
统 以 减少 不 断 变 化 的 大 气 辐射 ， 并 在 主 焦 点 安装 了 一 个 冷却 接收 器 ， 带 宽 为 3GHz， 中 心 频 
率 为 13GHz。 该 望远镜 对 每 个 射电 源 监测 频次 大 约 是 每 周 监 测 两 次 。 伽 马 射 线 数 据 来 源 于 
孔径 测 光 法 ( 测 光 半径 为 1” ， 最 大 天 项 角 105”)， 因 此 并 没有 去 掉 伽 马 射线 背景 光子 ， 但 
是 背景 光子 流量 比较 稳定 ， 对 我 们 的 光 变 分 析 的 影响 可 以 忽略 不 计 “〈 我 们 对 J0929+5013 和 
J2146-1525 伽 马 波段 的 数据 用 费 米 工 具 包 处 理 , 得 到 7 天 bin、10 天 bin 和 30 天 bin 的 数据 。 
用 LSP 方法 和 加 权 小 波 Z 变换 方法 计算 得 到 结果 基本 相同 ) 因 此 本 文选 取 30 天 bin 的 数据 。 
我 们 采用 光 变 指数 4 判断 天 体 的 活跃 程度 19]; 
A= J Smin . 
Sai Siin 
(1) 式 中 fiwax 表示 流量 的 最 大 值 、fiwin 表示 流量 的 最 小 值 。 通 过 计算 可 以 得 到 J0929+5013 和 
J2146-1525 射电 波段 的 光 变 指数 分 别 为 0.55 和 0.35， 伽 马 波段 的 光 变 指数 分 别 为 0.53 和 
0.46。 我 们 知道 光 变 指数 越 接 近 1， 天 体 越 活跃 ， 这 表明 J0929+5013 和 J2146-1525 射电 和 
伽 马 波 段 比较 活跃 。 图 1 是 J0929+5013 射电 和 伽 马 波段 的 光 变 曲线 ， 图 2 是 J2146-1525 射 
电 和 伽 马 波 段 的 光 变 曲线 。 横 坐标 为 简化 儒 略 日 (Modified Julian Day, MJD)， 纵 坐标 为 流量 
及 误差 〈 大 部 分 误差 很 小 ， 图 中 不 能 清晰 分 辨 ) 。 
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Fig.1 The Light curves of CGRaBS J0929+5013 in radio and gamma-ray wave bands. 
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图 2 CGRaBS J2146-1525 在 射电 和 伽 马 波段 的 光 变 曲线 。 
Fig.2 The light curves of CGRaBS J2146-1525 in radio and gamma-ray wave bands. 
2 LSP 周 期 分 析 


漆 变 体 的 周期 性 光 变 研究 , 对 理解 其 内 部 辐射 的 物理 机 制 和 辐射 区 的 几何 性 质 有 重要 意 
义 ，Lomb-Scargle periodogram (LSP) 方法 是 基于 侍 里 叶 变 换 ， 使 用 正弦 函数 模型 ， 可 以 对 
非 均匀 采样 的 时 间 序 列 周期 图 进行 相位 修正 , 能 够 在 一 定 范 围 内 对 非 均 匀 采 样 的 时 间 间 隔 引 


起 的 误 报 周期 进行 修正 。 可 以 比较 好 地 寻找 隐藏 在 噪声 中 的 准 周期 振荡 。 其 基本 公式 0 为 


> (x(1;)—x) cos @(t, —7) > (x(¢,)—x)sin ot, —7) 
j=l 十 j=l 


N N 
>》 cos’ w(t; — 7) > sin’ a(t, -7) 
jal j=l 


p,(@) = 5 


DRP, t, 为 时 间 序 列 的 平均 值 , N 为 数据 点 的 个 数 ,z 为 相应 的 时 间 的 相位 修正 : 
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N 
> sin? 2ot; 
tan(2@T ) = E ___ ， (3) 
>) cos? 2at,; 
j=l 


我 们 对 LSP 方 法 计算 的 结果 用 蒙特 卡 洛 模拟 进行 置信 度 评估 10， 我 们 对 LSP 做 对 数 寡 率 
拟 合 ， 得 到 J0929+5013 和 了 2146-1525 在 射电 和 伽 马 波段 的 寡 率 指数 BUG HA B =0.5、 


B=0.1; B=0.7、B =0.3。 发 现 暴 率 指 数值 偏 小 ,我们 认为 主要 的 原因 可 能 是 受 白 噪声 影 


响 ， 在 相当 一 部 分 情况 下 ， 军 率 指数 大 约 为 1 0 。 
图 3 和 图 4 分 别 是 J0929+5013 和 J2146-1525 在 射电 和 伽 马 波段 的 蒙特 卡 罗 置 信和 度 分 析 结 
果 ， 其 中 绿色 线 表示 准 周期 频率 图 ， 峰 值 就 是 准 周期 频率 ， 黑 色 、 红 色 、 蓝 色 线 分 别 表 示 蒙 


特 卡 罗 模 拟 的 2 a (95%) 2.60 (99%). 30 (99.7%) #15 BE HH AE 
图 3 中 的 (a)、(b) 分 别 表示 J0929+5013 在 射电 和 伽 马 波段 的 蒙特 卡 罗 置 信 度 分 析 结 果 ; 通 
= 过 计算 得 到 J0929+5013 在 射电 波段 有 一 个 1435 天 〈 约 3.93 年 ) 的 置信 度 >3 o 水平 的 准 周期 。 
一 J0929+5013 在 伽 马 波 段 有 五 个 置信 度 >2 Ga 的 峰 ， 其 中 112 天 、146 天 、203 约 为 53.4 天 《卫星 
: 的 进 动 周期 ) 2 倍 、3 倍 、4 倍 ， 我 们 推测 它们 是 飞船 进 动 引起 的 高 次 谐 频 ，192 约 为 半年 ， 也 


是 一 个 谐 频 。 我 们 认为 311 天 可 能 是 分 马 波段 一 个 很 弱 的 可 能 准 周 期 ， 有 待 进一步 确认 。 


= 如 图 4 中 的 (a)、(b) 所 示 ,，J2146-1525 在 射电 波段 有 一 个 1321 天 〈 约 3.62 年 ) 的 置信 度 >3 o 
© 水 平 的 准 周期 。J2146-1525 在 伽 马 波段 有 三 个 置信 度 >2 o 的 峰 , 其 中 108 天 约 为 53.4 天 的 2 倍 ， 
= 130 天 约 为 月 亮 周 期 27.3 天 的 4 倍 。 我 们 认为 一 个 1388 天 是 伽 马 波 段 的 一 个 很 弱 的 可 能 准 周 
N 期 ， 有 待 进一步 确认 。 
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(a) (b) 
图 3 CGRaBS J0929+5013 在 射电 (8B =1D) 和 伽 马 (8B =1) 波 段 准 周期 震荡 的 蒙特 卡 罗 置 信和 度 分 析 结 果 。 


Fig.3 The confidence analysis results of quasi periodic oscillation of CGRaBS J0929+5013 in radio( B =1)and gamma-ray( B =1)wave 


bands by Monte Carlo simulation. 
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(a) (b) 
414 CGRaBS J2146-1525 在 射电 (8B =1.3) 和 伽 马 (8B =1) 波 段 准 周期 震荡 的 蒙特 卡 罗 置 信和 度 分 析 结 果 。 


Fig.4 The confidence analysis results of quasi periodic oscillation of CGRaBS J2146-1525 in radio( B =1.5)and gamma-ray( B =1)wave 


bands by Monte Carlo simulation. 


3 加 权 小 波 Z 变 换 周 期 分 析 


小 波 分 析 (Wavelet Analysis) 是 一 种 由 傅 里 叶 变 换 演 化 而 来 的 分 析 方 法 ， 能 将 时 域 和 频 域 
局 部 化 ， 通 过 平移 和 伸缩 变换 可 以 把 小 波 母 函数 变 为 小 波 函 数 03，(4) 式 中 a 和 b 分 别 表示 
伸缩 尺度 和 平移 参数 : 


1 t—b 
Walt) = aa ae ， (4) 


Foster 定义 了 加 权 小 波 变换 (WWT), (5) 式 中 N e 为 有 效 数据 点 个 数 ， 矿 为 数据 的 加 权 偏 差 ， 


为 模型 函数 的 加 权 偏 差 ， 


eff y 
2y 


x 


_ (Ny -DV 


WWT (5) 


TH 


这 原 至 


在 低频 部 分 分 辩 率 很 高 ， 导 致 低频 部 分 有 效 数 据 Ns 多 于 高 频 部 分 ， 所 以 WWT 
的 值 倾向 于 高 频 部 分 ， 因 此 Foster 按照 Z 统计 量 对 加 权 小 波 Z 变换 n9 进 行 了 定义 : 
_ Ney 3) 六 


3 


2y -V,) | 


我 们 根据 加 权 小 波 Z 变换 的 基本 原理 , 编写 python 程序 处 理 J0929+5013 和 J2146-1525 
在 射电 和 伽 马 波段 的 数据 。 图 5 和 图 6 中 黑色 曲线 作为 周期 的 函数 图 ， 其 中 红色 、 绿 色 、 蓝 
色 线 分 别 表示 基于 红 品 声 的 2o 、2.60 、3o 置信 概率 线 , 黄色 虚线 表示 基于 白 噪声 的 99.7% 
置信 概率 线 。 根 据 加 权 小 波 Z 变换 图 中 的 极 大 值 可 以 确定 光 变 准 周 期 。 
图 5 中 (a)、(b) 分 别 表示 J0929+5013 在 射电 和 伽 马 波段 用 WWZ 方法 计算 的 结果 ，(a) 
图 显示 J0929+5013 在 射电 波段 存在 一 个 置信 和 度 >3o 的 极 大 值 1435.0 天 A 3.93 Æ), Ab) 
图 显示 J0929+5013 在 伽 马 波段 存在 一 个 置信 和 度 >2o 的 极 大 值 306.7 天 ( 约 0.84 年 ); 


(6) 


A 


J0929+5013 在 射电 和 伽 马 波段 基于 白 噪 声 的 置信 度 都 >3 o ， 这 是 一 个 比较 强 的 假设 。 我 们 
认为 J0929+5013 在 射电 波段 存在 一 个 1435.0 天 的 比较 可 信 的 准 周期 ， 在 伽 马 波段 306.7 天 
的 准 周期 还 需 进一步 确定 。 

如 图 6 中 的 (a)、(b) 所 示 , J2146-1525 在 射电 波段 存在 一 个 置信 度 >3 o 的 极 大 值 1320.5 RK 
3.61 年 ) ， 了 2146-1525 在 伽 马 波段 存在 一 个 极 大 值 1384.8 K (A) 3.80 年 ) ;了 2146-1525 在 
射电 和 伽 马 波段 基于 白 噪声 的 置信 度 都 >3 o 。 我 们 认为 卫 146-1525 在 射电 波段 存在 一 个 
1320.5 天 的 比较 可 信 的 准 周期 ， 在 伽 马 波段 1384.8 天 的 准 周期 还 需 进 一 步 确定 。 
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Fig.5 The WWZT test result for quasi periodic oscillation CGRaBS J0929+5013 at radio and gamma-ray wave bands. 
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Fig.6 The WWZT test result for quasi periodic oscillation search CGRaBS J2146-1525 at radio and gamma-ray wave bands. 


4 相关 性 分 析 


期 , 那么 同一 个 天 体 不 同 波段 的 光 变 有 没有 相关 性 呢 ? 为 了 研究 这 个 问题 , 我 们 用 离散 相关 
分 析 方 法 分 别 对 J0929+5013 和 了 J2146-1525 的 射电 和 伽 马 波 段 数 据 做 相关 性 分 析 。 离 散 相关 函 
数 (Discrete Correlation Function, DCF) 比 较 适 合 分 析 时 域 分 布 不 均匀 的 天 文 观测 数据 中， 基 
本 公式 为 : 


上 面 我 们 已 经 用 不 同 的 方法 计算 了 J0929+5013 和 J2146-1525 在 射电 和 伽 马 波段 的 准 周 


(a, —a)x(b,—b) 


O, 一 Op 


UDCF, = (7) 


(DRF a Mb 分 别 表示 两 组 数据 的 平均 值 ，o, 和 ca, 为 两 组 数据 的 标准 偏差 。 


DCF(t) = —_YUDCE, f (8) 


将 (7) 式 代入 (8) 式 中 ， 得 到 了 相关 系数 。 

图 7、 图 8 分 别 是 J0929+5013 射 电 和 伽 马 波段 、JP146-1525 射 电 和 伽 马 波段 的 相关 性 分 析 
结果 ， 从 图 中 可 以 看 出 相关 系数 DCF 的 最 大 值 仅 为 0.32 和 0.30， 相 关 性 在 0.5-1 之 间 我 们 才 认 
为 两 者 可 能 存在 相关 性 ， 因 此 J0929+5013 射 电 和 伽 马 波段 、J2146-1525 射 电 和 伽 马 波段 没有 


明显 的 相关 性 ， 我 们 认为 两 个 波段 的 辐射 可 能 来 自 不 同 的 区 域 
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图 8 CGRaBS JP146-1525 在 射电 和 伽 马 波段 相关 性 分 析 结 果 。 


Fig.8 The DCF analysis result of CGRaBS J2146-1525 in radio and gamma-ray wave bands. 


5 讨论 与 结论 


我 们 用 LSP 方 法 和 加 权 小 波 Z 变 换 方法 对 J0929+5013 和 卫 2146-1525 射 电 和 伽 马 波段 的 数 
据 进 行 分 析 ， 得 到 J0929+5013 在 射电 波段 存在 一 个 置信 度 >3 o 约 1435 天 的 准 周期 光 变 ， 在 
伽 马 波段 存 在 置信 度 >2 o 的 很 弱 的 309 天 左右 的 准 周期 ， 用 DCF 方 法 分 析 两 个 波段 的 相关 
性 ， 发 现 这 两 个 波段 没有 明显 的 相关 性 ， 我 们 认为 两 个 波段 的 辐射 来 自 不 同 的 区 域 。 
J0929+5013 在 射电 波段 存在 一 个 置信 度 >3 o 约 1321 天 的 准 周期 光 变 ， 在 伽 马 波 段 也 存在 一 
个 很 弱 的 1385 天 的 可 能 准 周期 光 变 , 同样 , 我 们 用 DCF 方 法 对 两 个 波段 的 数据 做 相关 性 分 析 ， 
分 析 结 果 表 明 这 两 个 波段 的 相关 性 也 很 弱 。 由 于 数据 点 有 限 ， 以 及 光 变 曲线 的 时 间 跨 度 不 够 
长 ， 大 约 只 有 所 测 QPO 时 标的 四 倍 ， 我 们 给 出 的 只 是 两 个 QPO 候 选 体 。 耀 变 体 中 QPO 是 非 
常 军 见 的 ， 因 此 需要 更 长 时 间 和 更 多 的 数据 来 证 实 。 

对 于 耀 变 体 的 准 周期 振荡 (QPO) 的 起 源 和 理论 解释 都 还 存在 争议 。 目前， 主要 的 解释 
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有 三 种 ;第 一 种 双 黑 洞 引起 的 螺旋 喷 流 模型 4081， 第 二 种 单 黑 洞 的 吸 积 盘 脉动 模型 49]， 第 三 
种 单 黑洞 引起 的 喷 流 进 动 模 型 名 。 我 们 基于 双 黑 洞 的 螺旋 喷 流 模型 对 J0929+5013 和 
J2146-1525 的 质量 做 了 估算 , 根据 Rieger 在 2004 年 2 给 出 了 观测 运动 准 周期 尸 与 实际 物理 


驱动 周期 P, 之 间 的 关系 。 
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~ Vs 
Pa * ie (8) 
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(8) sty, 为 体 洛 伦 兹 因子 约 等 于 7.52]， 书 的 单位 为 年 ， 我 们 用 比较 可 信 的 射电 波段 的 
观测 准 周期 ,它们 分 别 为 3.93 年 和 3.62 年 , 红 移 分 别 为 0.37 和 0.70。J0929+5013 All J2146-1525 


司 之 间 的 质 


的 实际 物理 驱动 周期 分 别 为 已 分 别 为 161.36yr 和 119.78yr WE WI 
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EHE R <13, 称 之 为 “ 主 并 合 ” (major merger) 。 如 果 3< R S104 ZA “RFE” (minor 
mergenP3。 主 黑洞 的 质量 估算 公式 为 P4， 
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as EDS 6 
M ~ PiR510 Mo (9) 
(9) 式 中 ,假设 主 次 黑洞 质量 比 为 R=3/2。 把 参数 代入 (9) 式 ,得 到 J0929+5013 和 J2146-1525 
主 黑洞 质量 MM 分 别 约 为 4.3x10” Ms M 2.7x10°M , 同 理 ， 如 果 使 用 y=15R5 作 为 参数 ， 


并 合 的 SMBBH 系统 的 主 黑洞 质量 分 别 为 4.0x101° Ms 和 2.5x1019 Moo WE J0929+5013 


和 J2146-1525 主 黑洞 质量 M 分 别 为 4.3x10?M 和 2.7x10” Mo 
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Correlation and Quasi—Periodic Oscillation Analysis of 
Light Variation in Blazar CGRaBS J0929+5013 and 
J2146-1525 


Chen Junping’’, Ma Li', Gong Yunlu', Chen Zhihui', Yang Xing”’, Chang Xin’, Li Yong’, Yi 
Tingfeng… 
(1.School of physics and electronic information, Yunnan Normal University, Kunming 650500, China, yitingfeng98@163.com; 2. 
Guangxi Key Laboratory for the Relativistic Astrophysics, Nanning 530004, China; 3.School of Physical Science and Engineering 


Technology, Guangxi University, Nanning 530004, China) 


Abstract: We used the Lomb-Scargle periodogram (LSP) and the weighted wavelet Z-transform 
(WWZ) methods to analyze the observational data of the blazars CGRaBS J0929+5013 and 
J2146-1525 at the 15 GHz radio and gamma-ray bands for about 12 years. The result of analysis 
showed that the blazars CGRaBS J0929+5013 and J2146-1525 probably have a quasi-periodic 
oscillations (QPOs) of about 1435 days (3.93 years) and 1321 days (3.62 years) in the radio band, 
respectively, with a confidence level of >3 o (99.7%). However, the periodicity of the light curve 
of the gamma band is weaker. Based on the spiral jet model in the supermassive dual black hole 


system, the primary black holes masses M of the blazars CGRaBS J0929+5013 and J2146-1525 


are estimated to be 4.3x10? M, and 2.7x10? M , , respectively. In order to explore the origin of 


optical variation in radio and gamma-ray bands, we calculated the correlation between the blazar 
CGRaBS J0929+5013 and J2146-1525 radio and the gamma band using the discrete correlation 
function, and found that there is no obvious correlation between them. This result indicate that the 


radiation may come from different areas. 


Keywords: Blazar, CGRaBS J0929+5013; CGRaBS J0835-6835; LSP Method; WWZ Method 


